Collective Rabi Oscillations and Cold Collisions. 
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We examine cold atomic collisions within a resonant optical cavity. The quantized cavity mode 
can be used to manipulate the collisions between the cold atoms, such that periodic exchange of 
excitations between the atoms and the electromagnetic field strongly alters the collision dynamics. 
A colliding pair of atoms can thereby oscillate between its ground and excited states during the 
collision time. Using a semiclassical model, it can be predicted that such Rabi-like oscillations are 
revealed in the atomic trap-loss probabilities, which show maxima and minima as a function of the 
detuning between the frequencies of the mode and the atomic transition. 



Exchange of excitations between the energy levels of 
atoms or molecules and the quantum electromagnetic 
field ||l| is one basic interaction process between matter 
and light. In this context, the Jaynes-Cummings model 
1^ , describing a single two-level system and a monochro- 
matic lossless field, reveals several characteristics of this 
interaction. In the quantum Rabi oscillation in 
particular, a single energy quantum is periodically ex- 
changed between system and field. On the other hand, 
when several atoms or molecules interact with the same 
field, quantum coherence can build among them leading 
to well known collective behaviors, such as superradiance 
|^|-^, for which multiparticle entanglement plays a major 
role. It should be noticed, however, that the systems with 
which the field interacts need not be solely composed of 
stable atoms or molecules. The above characteristics of 
the matter-field interaction can be present as well when 
two cold atoms collide in presence of a quantum field. 
An additional effect is that the very dynamics of the cold 
atomic collision ||l^,|ll| can be strongly modified by this 
field. Indeed, pairs of colliding atoms so far apart from 
each other that their direct mutual interaction is negligi- 
ble, can be entangled by the field and, thereby, influence 
one another in a nonlocal way. In the following we study 
cold atomic collisions within a gas of cold atoms trapped 
in the center of a high-Q cavity. It is possible to show 
that an analogous collective Rabi oscillation can show up 
in the trap-loss probabilities as a function of the detun- 
ing between the cavity mode frequency and the atomic 
resonance. 

The internuclear potential energy of two colliding 
atoms depends on the electronic states. In particular, the 
inverse-cube law dipole-dipole potential ±1/R^ between 
alkali neutral atoms, separated by a distance R, predom- 
inates when the asymptotic atomic states involved are 
n5i/2 and nPi/2 (Fig. p. When the colhsion is slow 
| pO| , pl| , the atoms can undergo changes in their electronic 
states during the collision, either as spontaneous decay 
or induced transitions. In a previous work |jl^, it was ex- 
amined how cold collisions can be manipulated as spon- 
taneous decay is driven by the colored vacuum of a cavity. 
The highly increased emission rate of multiply entangled 
pairs being able to emit coherently to the same cavity 



mode was then predicted to practically interrupt the col- 
lision process. In the present case, we neglect the cavity 
loss (high-Q) and allow reabsortions of the field energy. 
We study how these induced transitions affect cold col- 
lisions. Cavity Quantum Electrodynamics (CQED) ef- 
fects on atomic motion in high-Q cavities have been in- 
vestigated recently. Modifications of mechanical forces 
of light acting on atoms [|l3|,Q, or appearance of new 
quasi-bound molecular states of two colliding atoms |l^] 
illustrate the interplay between CQED and cold atoms. 




Internuclear separation (R) 
FIG. 1. Excited state long range dipole-dipole attractive 
potential U = —Cs/R'^ and the ground state van der Waals 
potential 1/R^, as functions of the internuclear distance R. 
The energy difference of the latter to the asymptote of U is 
the atomic separation Hlja between nSi/2 and nPi/2- 

A pair of colliding atoms under the attractive poten- 
tial UiR) — —Cs/R^, with C3 a constant, is described 
as in by a two-level system. Its energy splitting 
hujji — hujA + U{R) depends on the internuclear distance 
R and approaches the atomic energy difference TujJa be- 
tween nS 1/2 and nPi/2 as R — > 00. In this asymptotic 
limit, the excited state of the pair denoted by |e) cor- 
relates to nSi/2 + nPi/2, whereas the ground state de- 
noted by \g) correlates to nS'1/2 + nS 1/2 and the van der 
Waals potential in this latter case is neglected. Since 
the atoms are weakly bound by U (i?) in |e) we call them 
quasimolecules, extrapolating this denomination to \g) as 
well. The formation of this quasimolecule occurs at the 
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Condon point Rc, the distance at which lor becomes 
resonant with an excitation laser of frequency lul < lua- 
The nuclei may then attract each other starting from 
Rc- An eventual spontaneous decay from |e) to \g) re- 
leases a kinetic energy AK — Tilol — TT-oJ-y (see Fig |^). 
If AK is more than twice the trap depth Vbj i-e-, if the 
decay occurs below Re, the pair is ejected from the trap 
in the so-called radiative escape collision This 
two-level approximation neglects hyperfine splittings of 
the atomic fine structure and possible (anti) crossings of 
potential curves, whereas the nuclear motion is treated 
semiclassically. Nevertheless, this turns out to be reli- 
able for a satisfactory description of cold collisions in the 
case of ®^Rb, a range of hundreds of MHz for the detun- 
ing 5 = ujl ~ i^A, and a weak excitation laser intensity 
pO| , pl]Jl6[[r^ . Trapping the atoms in their ground state, 
e.g. with a far off-resonance trap (FORT) ||l8|, is most 
suited since optical pumping effects by the trapping light 
can be avoided; ejected atoms can then be controlled 
more efficiently by a separate excitation light field. 

For a high-Q cavity, the main difference with our pre- 
vious work is the exchange of excitations between the 
cavity mode and the quasimolecules. A set of N identical 
quasimolecules is described by 



If the excitation laser beam is normal to the cavity 
axis, the quasimolecules are excited independently from 
each other since the wavelength 2Tic/ix}R is in the opti- 
cal domain and is much shorter than the average sep- 
aration between quasimolecules | [l9[ |. This corresponds, 
for each quasimolecule, to a no-cavity condition since, 
in the optical domain and for alkalis, typically <C F, 
where F is the spontaneous decay rate of state |e) [||. 
In contrast, using for the excitation the quantized mode 
of the cavity by injecting the laser through the cavity 
axis (with ujl = ^c), the quasimolecules become in- 
distinguishable to the mode a^a and all of them have 
a quantum probability amplitude to be excited. They 
therefore end up in a multiparticle entangled excited 
state. Considering only one excitation, it follows from 
Hint that this state is 0) = /Cl, where 

|1,0) = \eg---gm, |2,0) = l^e • • • .9, ) |0) and so forth, 
with energy Ee{R) = hajR. It couples to the ground 
state |G, 1) = \g ■ ■ ■ g)\l), with energy Eq = HlJc of one 
field excitation in the mode and all quasimolecules in 
their ground state \g), by a coupling constant 
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where each Uzi is a Pauli spin matrix. Due to the nuclear 
motion, lor becomes formally time dependent since R 
may change with time. The actual value of N depends 
on the detuning 5, the total number Na of atoms and 
their density ua- The hamiltonian of the quantized field 
of the cavity mode is in turn given by 



(2) 



where ujc is the cavity resonance frequency and (a), 
the bosonic creation (annihilation) operator of the field. 
The interaction hamiltonian in the rotating-wave approx- 
imation is 
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These two states form a closed subspace for the total 
hamiltonian 



H — H„i + He + Hi- 



(6) 



It may be asked whether the evolution of the states 
\E, 0) and |G, 1) under this hamiltonian starts imme- 
diately after one excitation is injected into the cavity 
mode. Consider the case of one single pair of free atoms 
approaching each other in the state \g) (asymptotic to 
5'i/2 + 'S'i/2)- In the presence of a light field, the pair may 
adiabatically pass from \g) to |e). If a Landau-Zener ap- 
proximation [E^,^ holds, the probability for this to oc- 
cur is 1 — e~^^^, where A = h^l'^ /voo\U'{Rc)\, ^ being 
the coupling (Rabi frequency) between \g) and |e), and 
Voo , the asymptotic relative velocity. This approximation 
is satisfactory for detunings \5\ ^ lOF/i and small inten- 
sity ll^. Analogously, the probability Pe to excite 
the collective state 0) from |G, 1) may be estimated 

by 



Here, the Pauli matrices crj and Ui are raising and low- 
ering operators, respectively, acting on |e) and \g) of the 
i-th quasimolecule. The individual Rabi frequencies 



Vli = 2£{uje) fci^i) e • di/h 
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depend on the field strength per photon £{uj) ~ 
{V being the mode volume), polarization 
e, mode profile /c(r), and the molecular dipole moment 
di of the transition \g) — > |e) (whose absolute value 
|di| = -\/2|dA|, being the atomic dipole moment of 
nSi/2 nPi/2). 
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where the collective coupling fl replaces 51. Since f2 > 17, 
an immediate consequence of the cavity mode is to in- 
crease the excitation probability. Once |£',0) is formed, 
the quasimolecules' evolution is then governed by Eq.(^. 
Actually, we use a large value for Q and the Landau- 
Zener approximation may fail Nevertheless, with 
equal fractions of excited quasimolecules in both situa- 
tions with and without cavity, the excitation probability 
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can be factored out. Therefore, only the trap loss prob- 
ability that an excited quasimolecule, either in the col- 
lective state or in |e), is ejected from the trap becomes 
significant. 

The most notable feature arises as time goes on. One 
expects a Rabi oscillation between \E, 0) and \G, 1) as the 
excitation is exchanged between the quasimolecules and 
the cavity mode. This implies that the quasimolecules 
switch between |e) and \g) at the collective rate Cl, al- 
tering completely the collisional encounter compared to 
how it would proceed without the quantum field; and the 
larger the number of quasimolecules N, the faster such 
oscillation will be , since fl scales with ViV. This oscil- 
lation looses strength as ujr detunes gradually from Uc- 
When \ujji — ujc\ ^ ^, we may neglect Hint and let the 
quasimolecules evolve freely, subject at most to sponta- 
neous decay. It is interesting to mention here that the 
resonance condition is determined by Q instead of the 
cavity (or molecular) linewidth. This is analogous to a 
power broadening effect, even though we consider a single 
excitation in the cavity: the "power" here is increased via 
an increase of the collective dipole. If there is a proba- 
bility that the quasimolecules remain excited on entering 
the off-resonance region R < R' , determined by 

(^Rc -i^R' (8) 

they may reach R < Re where spontaneous decay of 
\E,0) ejects a pair from the trap. In contrast, if the in- 
teraction time is such that the quasimolecules leave the 
resonant region R' < R < Rc in the ground state |G, 1) 
(having thus undergone a so-called 7r-pulse), practically 
no trap loss will be recorded at all. From the dependence 
of this interaction time on the detuning S, we deduce in 
the following an expression for the trap-loss probability 
as a means to observe and measure this collective Rabi 
oscillation of colliding cold atoms as a function of 5. 

If tc ^ 27r/ri is the interaction time spent in the res- 
onance region R' < R < Rc, we keep the resonance 
condition ujn « ujc, since \ujr — uJc\tc ^ x 27r/f2 ~ 2tt 
completes at most a single cycle. Analogously, the time 
variation of ujr ^ is neglected: the phase 6{tc) = 
Jg" dr \cJii\T that such time dependence would introduce 
to the dynamics is at most of order 27r if we use that 
tc ^ 2TT/il and lAwflJ ^ (l. A full solution is possible, 
but this suffices for our purposes, limiting us, however, to 
few oscillations. The spontaneous decay rate of the state 
|e) is r = 2Ta, where I/Ta is the lifetime of the atomic 
state nPi/2 and retardation effects can be neglected 
{R is safely smaller than 2ttc/ujr). Since the cavity mode 
solid angle is small, the rate of emission into the rest of 
free space is nearly F. For the quasimolecules and the 
cavity mode, the Liouville-von Neumann equation with 
dissipation is 

1 ^ F 



where the cavity dissipation is neglected. The dissipa- 
tor in this equation assumes an incoherent decay of each 
quasimolecule independently from each other, since the 
emitted wavelength is shorter than the average separa- 
tion of quasimolecules When the emission is to- 
wards the cavity mode, in contrast, a collective effect in 
the quasimolecules' emission is fully contained in the in- 
teraction with the field. 

For the initial condition p(0) ~ \E,0){E,0\, we calcu- 
late the probability 

PE{t) = tv[\E,0){E,0\p{t)] (10) 

that the system remains in the state \E, 0) after a time 
interval t. The solution can be found by enlarging the 
subspace spanned by \E, 0) and \G, 1) to include the vac- 
uum state \ V) = jG, 0) = \gg ■■■ g)\0) to which the system 
evolves due to dissipation. In the interaction representa- 
tion, one has 

= Pe\E, 0) (£;, 0| + pg\G, 1) (G, 1 1 + pv\G, 0) (G, 0| 
+ {ceg\E,0)(G,1\ + cev\E, 0)(G,0| 

+CGy|G,l)(G,0|-l-adj.) . (11) 

The time-dependent probabilities pe, Pg, and pv, and 
the coherences ceg, cev, and cgv are found by sub- 
stitution of pi(t) into the interaction representation of 
the Liouville-von Neumann equation and projecting into 
each component of pi(t). The damping terms generate 



no new state and we get for the p's 

PE = -^PE + i^{cEG - C*eg), (12) 
PG = -i^icEG - C*eg), (13) 
PV^^PE, (14) 

and for the c's 

F 

CEG^~^CEG+i^{PE~PG), (15) 

F 

CEV ^ -^CEV - i^CGV, (16) 

CGv = -i^CEV- (17) 



From the last two equations, it follows that CEv(t) — 
CGvit) = if this holds initially. From Eq.(|l|), Eq.(yl), 
and Eq.(p^), pE{t) can be solved straightforwardly. Since 
> F/4 (condition for Rabi oscillations) is supposed to 
be fulfilled with the parameters we adopt, the solution is 
the underdamped one 



Pn{t) =PE{t) = 




where 

f}= - (F/4)2, n>r/A. (19) 
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As already mentioned, this solution is to be applied so 
long as the resonance condition — Wc| ^ ^ holds, 
that is, as the quasimolecules-field interaction is signifi- 
cant. For jcjfl' — ujc\ ^ ri, on the other hand, we neglect 
this interaction. The solution then is formally equivalent 
to taking the limit Q — > in the overdamped regime, 
which follows from Eq.(|lJ) by substituting — i|/3| for /? 
(i.e., solving Eq.(§) without Hint), 



PE{t) 



(20) 



As expected, this equation implies the value F for the 
decay rate of the collective state \E,Q), the same as the 
molecular value for |e). This can be derived directly from 
a perturbation calculation a la Fermi's golden rule by 
coupling the quasimolecules to all the electromagnetic 
vacuum modes except the cavity mode Uc, whose small 
solid-angle is neglected. It reflects the fact that the aver- 
age separation between quasimolecules is larger than the 
wavelength of the emitted radiation, and the indistin- 
guishability of which quasimolecule emits into free space 
can no longer hold |l9|] . 

We need pn to obtain the probability that one pair 
of atoms is ejected from the trap after excitation by the 
cavity mode. The transition from \E,0) to |G, 1) should 
then occur for R < R^. In the first passage in this region, 
the probability li that a pair of atoms from any one of 
the N quasimolecules is ejected is constructed as 



N 



Po(ic)e-*'^(l 



-2t,T 



(21) 



where {f , or tc) is the time interval spent between 
R = and Re {Re and R' , or R' and Rc)- This h is thus 
a sum of products of conditional probabilities jlOj , which 
in our case involve Eq. ([l8|) and Eq. (po|) . Now, by letting 
R evolve further, the quasimolecules can vibrate between 
Rc and R — before emission takes place, that is, they 
may pass several times across Re while still in 0). 
Summing over these multiple passages [ p^ and compos- 
ing them with conditional probabilities, the probability 
Cc for a pair of atoms to be ejected at any time is then 



Cc = ll+l2+l3 + --- 
W ^ 2 



-tT 



n 

■ ± 

pn{tc) 
•••} 



(1 - e-'*^' ) 



-2(t'+te)r 



,-2(t'+i,)r 



Pn{2te) 
pn{2te) 



-,-t'r 



-tT 



sinh (F^e 



[e(t'+te)r_pj^(2iJe-(*'+*=)r] 



(l-e--=M 



(22) 



As in Eq. (|21|) , the summation over the N quasimolecules 
is factored since the multiple passages probabilities In 



are the same for each quasimolecule. In Z„, pn{tc) is the 
probability that a quasimolecule exits the region R' < 
R < Rc in the excited state |e) after having been excited 
at Rc, whereas po(2tc) gives the probability of crossing 
i?' in |e) from the right after having crossed it from the 
left in |e). 

The dependence of Cc on the collective Rabi frequency 
can be observed formally by changing the interaction 
time tc and thus making oscillate. This time interval 
in turn is a function of the detuning 6 = lUc — and 
can be obtained by integration of the energy conservation 
condition |lO| 



+ U{R) = const. . 



(23) 



Indeed, choosing f2 such that t' m (i.e., R' « Re, see 
Eq.(||)), we have TiVt = 2Vb (emission on resonance with 
the cavity will not lead to trap-loss) and neglecting the 
initial velocity R at R = Rc, it follows from Eq. (p3|) 



te^toiS) [l-f{S)], 



0, R' » Re) 



(24) 



where to = tc + tc is the total time interval between Rc 
and R = 0, and f{S) is the fraction of to spent between 
Rc and Re, 
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50 Jr 



du- 
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- 1 




(25) 

-1/3 



50 = 0.746 normalizing /(O) to unity. A slight compli- 
cation is the dependence of on (5 as well via the total 
number N = N{5) of quasimolecules, so that VL changes 
as tc{5) is changed. In order to obtain Q. — ^{5), we 
approximate for large N the summation in Eq.(||) by an 
average (e; being the orientation of di, and d its modu- 
lus) 

^■' = Y,m'' 

i 



= N- 



^^(l/c(r,:)P>(k-e,f) 
n 

f^d"^ 



(26) 



with |/c(ri)| ~ cos{zikc) {z along the cavity axis and 
kc = ijJc/c) and being randomly oriented, and denot- 
ing an (averaged) single quasimolecule Rabi frequency. 
The total number N can be estimated by counting all 
pairs of atoms whose separation R is such that ujr — uJc, 
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with a spread AR about the Condon Point Rc deter- 
mined by the linewidth T of the state |e), namely, 



N ■ 



Naua 4:TtR},AR 



1 
2 

1 , r 

-Naua '^t^Rcu 77^1 — 

2 \dujR/dR\Rc 



NaUa 



2ttC3 



(27) 



having used that C^/R^ = h\5\. 

For the sake of comparison, we calculate the trap loss 
probability Co that, in the absence of the cavity, |e) de- 
cays in the region R < Re after multiple passages across 
Re- Both Cc and Co are similar, except by pn{t) which 
is substituted for the pure decay pE{t) ~ e" 



-tr 



Co^ 



sinh (teT) 
sinh {tc + te)T ' 



(28) 



describing the trap loss probability of a statistical mix- 
ture of pairs of atoms colliding independently [|ll],f6l . 
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FIG. 2. Trap loss probabilities Cc and Co as functions of 
the detuning S — Uc — u>a (see text for details). Note the 
"oscillation" in Cc brought about by the interaction with the 
cavity quantized mode. The curve levels off and approaches 
Co for (5/27r < —800 MHz, where the interaction time tc goes 
to zero. 

For numerical estimates, we consider ^^Rb atoms and 
the atomic transition 5Si/2 5Pi/2- The wavelength of 
the transition is = l-KcjiOA = 795 nm, with atomic 
decay rate Fyi/27r = 6 MHz. The coefficient of the dipole- 
dipole potential is taken as C3 « 11 x 10~^^ erg [ p3|j24| 
and the detunning 5 = Wc — ^a < is made to decrease, 
starting from (5/27r = —350 MHz. For this (5, the Condon 
point is Rc — 366 A and the total time to = 1.07x 10~® s. 
The trap depth Vb is chosen of order 5 mK (~ 100 MHz) 
and would have to decrease as b is decreased in order to 
match the condition (see Eq.(ra) et seq.) 



This choice for allows for the quasimolecules to en- 
ter directly the trap escape region R < R^ as they exit 
the resonant region Re < R < Rc in which their in- 
teraction with the cavity mode is strongest. With such 
parameters, and S going down to ~ —1.0 GHz (when 
il/27T Ri 70 MHz), we can approximate /3 « and ne- 
glect the sine in Eq.dlT 



Pn{tc 



-t,r/2 



{nte) 



(30) 



For the detuning 5 = —350 MHz, one obtains tc = to/ = 
0.55to and te = {1 — /)to = 0.45to. The largest phase be- 
comes then nte w 27r X 200MHz x 0.55 x 1.07 x 10~^s = 
2.357r and Q tc will decrease as 6 decreases. In Fig.(|^), we 
plot Co and Cc, assuming the quasimolecules are equally 
excited in both cases so that the excitation probabili- 
ties can be factored out. Outside the range of detunings 
shown, the approximations in our model are no longer 
valid |l^jr^: for smaller |(5|, a full quantum mechanical 
treatement of the nuclear dynamics is required, whereas 
the discrete quantum vibrational levels of the dipole- 
dipole potential U{R) become resolved for \6\ ^ 1.0 GHz. 

Within this detuning window, nevertheless, it can be 
clearly seen the drastic change brought about by the cav- 
ity mode. The collective Rabi oscillations begin mildly at 
a large \6\ and become more pronounced at smaller val- 
ues of absolute detuning, where the interaction time tc 
is longer. At the minima of Cc, the quasimolecules leave 
the resonant region in the state \G, 1) (having undergone 
a TT-pulse) and the kinetic energy is insufficient to eject a 
pair of atoms from the trap. The maxima of Cc, where 
the quasimolecules leave the resonant region in the ex- 
cited state, are above Co since, without the cavity effect, 
the probability of a quasimolecule to reach the escape 
region in the excited state |e) is smaller. These maxima 
increase with decreasing 5 because the potential becomes 
steeper and, therefore, the time tc decreases. The effect 
of damping during the interaction with the cavity mode 
becomes less and less important and becomes closer 
to one (see Eqs.(|l|) and (|22|)). 

It should be noted that an actual realization of such 
an experiment is quite difficult. For a centimeter sized 
optical resonator with mirrors separated by Z « 1 cm, the 
mode volume V w irwll « 4.0 x 10~^ cm^, with a waist 
Wo — \J cljuic ~ 36 /im, so that using d = \/2 dA and 
F^ — Ad\Lo\/'ihc^ {loa ~ Wc), the averaged single quasi- 



2Vo = hVl, or fi/27r < 200 MHz . 



(29) 



molecule Rabi frequency in Eq. (|26|) is il/27r w 0.42 MHz. 
From Eq.([29|) and Eq. (p6|) , one needs at least N ^ 
(fi/ri)^ « 2.3 X 10^ quasimolecules (but less for larger 
\5\). Using Eq.(p7|), this could only be achieved with a 
gas of Na ^ 2.0 X 10^ atoms occupying all the mode vol- 
ume V , i.e., at a high density ua ~ 4.0 x 10^'^ cm~'^. A 
possible remedy that could do with less atoms at lower 
densities would imply, however, a smaller f2 and thus a 
substantially longer interaction time tc in order for the 
phase f2tc to reach ^ 27r and show a "coUisional Rabi os- 
cillation" ; for these longer times, the detuning would fall 
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below the limit \6\ ^ 10 Ta and our semiclassical approx- 
imation for the nuclear motion would break down |17| . 

In conclusion, it is a remarkable effect that the outcome 
of a collision encounter between two cold atoms can be 
so deeply affected when it takes place in the presence of 
a quantized cavity mode. The slowness of the collision 
makes possible the exchange of excitation between the 
colliding pair of atoms and the cavity field. Depending 
on how long the interaction is effective, the potential en- 
ergy of the collision can simply be "turned off" as the 
quasimolecule leaves the excitation within the cavity; in 
contrast, the quasimolecule may be taken closer to the 
radiative escape condition if this same excitation is re- 
turned back to it. In this process, several quasimolecules 
get quantum mechanically entangled to each other. The 
trap-loss probability as a function of the detuning can 
thus show a Rabi-like collective oscillation. 
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